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Simple algorithm of dynamics of Ising magnetic is described. The algorithm can 
be implemented on conventional digital computer and can be used for construction 
of specialized processor for simulation of ferromagnetic systems. The algorithm gives 
a simple way to calculate ID correlation functions for ID Ising magnetic. 



I. 



INTRODUCTION 



In recent past the analog computers were the only tool for solving the mathematical 
problems which could not be solved analytically. The analog computer is a device (typically 
- an electric circuit) whose temporal dynamics is described by equations similar to those 
which are to be solved. The analog computers have the virtue of being simple. The main 
defect of analog computers is absence of universality. For this reason for solving of any 
particular problem one must design a special analog computer suited for solving only this 
particular problem. The impressive achievements of semiconductor technology of recent 
decades made it possible to create a digital universal computer based on microprocessor - a 
programmable electronic device which is millions times more complicated than any of analog 
computers but allows one to solve wide range of problems just by entering of appropriate 
program. Despite the great power of modern digital universal microprocessors one can 
point out some problems whose solution require excessive time or even impossible. Having 
in mind one of these problems one can put the following question: is it possible (using 
the fantastic facilities of up-to-date technology) to create a specialized device (specialized 
processor) suited for solving only this particular problem? Imagine that we have constructed 
the specialized processor for solving of many-particle problem of atomic physics. In our 
opinion despite the loss of universality this processor would be of great interest. Below we 
suggest some algorithm for solving the problem of phase transition in Ising magnetic. This 
algorithm can be directly implemented by means of conventional digital computer. On the 
other hand in our opinion it is quite possible to create the specialized processor working 
according to this algorithm but much faster than conventional universal processor. This 
processor being much more simple device than the universal processor allows one to solve 
the problem of phase transition in Ising system and to obtain an arbitrary values of interest 
(energy, magnetisation, heat capacity). We estimate the complexity of this processor to be 
comparable with that of modern memory devices. 

Let us consider a trigger - scheme with two stable states which we denote by ±1. This 
trigger and spin 1/2 have much in common and below we will not differ these two. Suppose 
this trigger can change its state only when clock pulse coming. We consider the train of 
clock pulses to be equidistant in time. Consider the probability for trigger to change its 
state (the probability of corresponding spin flip) be depending on the trigger's state before 
clock pulse coming. Consequently the probability p+ for trigger to switch from state s = +1 
differs from that p- to switch from state s = —1: p+ ^ p^. Let these probabilities obey: 




s = 



±1. 



(1) 



Equation (1) has the sense of detail equilibrium principle for spin 1/2 in magnetic field h 
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and in contact with thermostat with inverse temperature (3. Since p± < 1 the normahzation 
constant ae must obey: 

ae < exp (2) 

If we now consider the statistical ensemble of such triggers then the kinetic equations for 
average number n± of triggers in states ±1 have the form: 

An; = -n%p+ + nip^ (3) 

An!_ = -n'_p- + n\p+ 

where An!,. - is increment of number of triggers in state s = ±1 after i-th clock pulse 
coming, n]^ - number of triggers in state s = ±1 to the moment of i-th clock pulse coming. 
The steady state solution of this equation has the form: 

1^ = ^ = exp(-2/3/.) (4) 

This corresponds to thermal equilibrium state of spin system in the external magnetic 
field h. 

Now let us consider the case of Ising magnetic i.e. the lattice comprised of N spins 
(triggers) coupled to each other in such a way that the magnetic field acting on the arbitrary 
spin is defined by configuration of the rest spins in the lattice. In the simplest case the 
magnetic field hi acting on spin with number i is produced by its nearest neighbours and if 
we denote the set of nearest neighbours by nn{i) then: 

hi^W J2 Sr, (5) 

r€nn(i) 

here the constant W characterize interspin coupling. By the analogy with the aforesaid let 
us consider the following dynamics of this system. The clock pulses act sequentially on all 
the spins (triggers) in the lattice - we call this round trip. During the round trip any spin 
(say i-th) may be overturned with probability defined by equation (1) with magnetic field hi 
defined by formula (5). We are interesting in the dynamics of this system under the action 
of the train of round trips. For the above dynamics of this "Ising magnetic" (consisting 
of triggers controlled by clock pulses) we now will obtain the kinetic equation for density 
matrix and will show that its steady state solution corresponds to thermal equilibrium. By 
the analogy with real magnetic let us describe the state of our "magnetic" (consisting of 
triggers) by the wavefunction whose i-th argument describe the state of i-th spin (trigger) : 

* = S2, ...Sat) (6) 
Introduce Oi - operator of i-th spin flip: 

Oi\Si, S2, Si, Stv) = S2, —Si, Sn) (7) 

Let us introduce the statistical ensemble of Ising magnetics and let o-(\&) be the number 
of magnetics in state \E' in this ensemble. Up to normalization factor the quantities (t(\1') 
are represent the diagonal elements of the density matrix of Ising magnetic. Let us consider 
Act' - the increment of a when clock pulse act on i-th spin (trigger) 
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Aa'i^) = - 5^ + ^ a{^)V^^^ (8) 

Here - is the probability of transition from ^'-state to $-state when clock pulse act 

on i-th spin. In accordance with the above dynamics of Ising magnetic the only non-zero 
probabilities are: 

(probability of transition from to $ = 0^^') = aeexp (j3hi{'^)si{'^)j (9) 
(probability of transition from $ = Oj^f to = aeexp ^/3/ij($)sj($)^ = (10) 

= aeexp f-/3/ii(^')si(^') 



The last equality follows from the fact that the field acting on i-th spin in states $ = Oj\& and 
\E' is the same while the value of i-th spin has the opposite sign (i.e. Sj($ = Oj^) — — Sj(^)). 
Then using equation (8) one can see that 



(11) 



a(Oi^) exp - /3/i,(*)si(*) j - a(^) exp \^(3hi(^)si{^ 
Let us show that this equation has the steady state solution in the form: 

ae,(*) = exp (Ai/(*)) (12) 



where H - is the Hamiltonian of Ising magnetic: 



^=YE E (13) 

j=l a=nn(i) 



We need to calculate creg($ = Oj^'). To do this note that: 



-ff(^)-2/i,(^)s,(^) (14) 



a=nn{i) * 

Now calculating (Teg($ = Oi'if) by equation (12) and substituting the result in to equation 
(11) it is easy to see that when A = — /3 the right part of equation (11) vanishes. So we 
see that (12) is the steady state solution of (11) and represent thermal equilibrium density 
matrix of Ising magnetic. Thus the above algorithm of sequential round trips prepare the 
system of coupled triggers in thermal equilibrium state. For transition probabilities be less 
than unit the value of ae should obey: 

<Si<exp(^- pm\W\y (15) 

here m - is the number of nearest neighbours. In our opinion it is possible to create a 
specialized processor working in accordance with this algorithm. The described algorithm 
can be implemented on conventional digital computer. In this should take the 

maximum possible value ae = exp ^ — /3m|VF|^ for system to relax as fast as possible. 
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In relaxed system one can observe magnetisation 5* = J^i^i, energy (13), heat capacity 
c = d{H)/dT, an arbitrary correlation functions. 

To demonstrate the aforesaid algorithm we simulate the phase transition in two- 
dimensional Ising magnetic by means of conventional computer. Fig.l shows the temperature 
dependence of heat capacity (top), magnetisation (middle) and energy (bottom) calculated 
for the case of = — 1 (ferromagnetic). Calculations were performed for lattice with sizes 
1000 X 1000. The procedure was as follows. At the beginning the system was prepared 
in the state with magnetisation close to its ultimate value S = N and with temperature 
much lower than the temperature of phase transition Tg. After that the above algorithm 
started with gradually increasing temperature 1//5. When temperature becomes close to Tc 
the magnetisation vanishes and heat capacity takes its maximum value. The value of Tc 
obtained in our calculations is in agreement with the exact formula of Kramers and Wannier 

i. 

The similar calculations can be performed for the case of zero initial magnetisation S = 
0. In this case total magnetisation is zero for all temperatures. Fig. 2 shows the spatial 
distribution of 2D magnetisation below (top picture, the domains are clearly seen) and 
above Tc (bottom). The heat capacity and energy temperature behaviour is similar to that 
in fig.l. 

The algorithm described can be directly generalized for the case of Ising system with an 
arbitrary interspin interaction W{r): 

ik 

To do this one should use the effective magnetic field in the form 

r 

instead equation (5). 

Possibly the described algorithm may be useful for checking the gauge theories of critical 
phenomena [1, 0, In this case the duality of some gauge models with respect to Ising 
system is exploited. The similar algorithms were described in |1, 0, HI. 

II. ID- CORRELATION FUNCTIONS. 

To check the above algorithm let us consider the exactly solvable one- dimensional Ising 
magnetic with nearest neighbours interaction. This problem was solved by Ising p| but above 
algorithm provide a simple way to obtain formulas (27) for correlation functions which are 
not very popular. In the case of ID Ising magnetic with nearest neighbours interaction the 
field acting on the i-th spin can be calculated by formula (5) as: 

hi = W{si_i + Si+i) (16) 

Suppose the clock pulse acts on i-th spin (trigger). Let us calculate the increment A(sj/), 
where / is an arbitrary function of of all spin variables except Sj. Multiplying both parts of 
equation (11) by Sif and summing over all states \l/ we obtain: 



A{sJ) = ae -J2(T{'^)fsiexp[(3W{si+i + s,_i)s,] 



(17) 



Passing from summation over ^ to summation over in the second sum (sj should be 
replaced by — Sj), denoting 

(3W = e, (18) 

and using the relation 

expas = cha + s sha, s = ±1, (19) 

we obtain: 

A{sJ) = ~2^{fs, cxp[/3iy(si+i + Si.r)si\) = (20) 
= -2ae(ch2^(s,/) + ^sh2^((si+i/) + (s,_i/)) + sh^e{si+,SiSi-if) 
Hence in the equilibrium state: 

ch^e{s,f) + ^sh2^((si+i/) + (Si_i/)) + sh^e{si+,SiSi-,f) = (21) 
Now we use this relationship to calculate the equilibrium correlation function: 

hp = . (22) 

This function depends only on the difference p of its indexes. Let / = Sj+p in equation 
(21). Then we have: 

A;pch^6' + ^ (fcp+i + A;p_i)sh26' + sh^6'(sj+iSjSj_iSj+p) = (23) 

To calculate the correlation function entering the last term let / = in equation 

(21). We have: 

ch.'^ e{si+iSiSi_iSi+p) + i {kp+i + kp_i)s\i2e + kps\i^e = (24) 

Hence: 

(sj+iSiSj-iSj+p) = -kptli^O - [kp+i + kp^i\th.O. (25) 
By substituting (25) in to (23) one can obtain: 

A;p = e(Vi + Vi), e = -^th(2^) (26) 

The solution of this equation under condition ko — l has the form: 

kp — exp(Q!p), C > 0) (ferromagnetic) (27) 

kp — (— l)^exp(Q;p), C < 0) (anti-ferromagnetic) 

with a (for p > 0) being the negative root of the equation: 

cha = ^ (28) 
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Formulas (27) were verified by direct computer simulation according to the above algo- 
rithm for ID Ising magnetic. 
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